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ABSTRACT 

Optical images of SN 1987A show a triple ring structure. The inner (dust) ring has recently increased 
in brightness and in the number of hot spots suggesting that the supernova shock wave has collided with 
the dense pre-existing circumstellar medium, a scenario supported by radio and X-ray observations. 
Such a shocked environment is widely expected to result in the acceleration of charged particles, 
and the accompanying emission of very high energy gamma-rays. Here, we report the results of 
observations made in 2004 and 2006 which yield upper limits on the TeV gamma-ray flux, which are 
compared with a theoretical prediction. In addition, we set upper limits on the TeV flux for four 
high energy objects which are located within the same field of view of the observation: the super- 
bubble 30 Dor C, the Crab-like pulsar PSR B0540— 69, the X-ray binary LMC X-l, and the supernova 
remnant N157B. 

Subject headings: gamma rays: search — supernova: individual (SN 1987A) 



1. INTRODUCTION 

SN1987A was the closest observed supernova in 380 
years, and the evolution of the remnant of this super- 
nova has been studied in great detail over the past twenty 
years. The radio intensity is growing and its rate of in- 
crease is increasi ng, with the spectral index being ob- 
served to flatten (|Stavelev-Smith et alj [20071. Imaging 
at radio wavelengths has enabled the expansion of the 
inner ring to be traced, and is revealing increasing struc- 
ture in the inner ring. The X-ray fluxes observed with 
XMM an d Chandra also continue to rise alm ost expo- 
nentially (jHabert et alj|2006t IPark et al.ll2007fh Around 
^4000 days after the supernova the X-ray light curve 
increased dramatically, attributed to the arrival of the 
supernova blast wave at the equatorial ring of circum- 
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stellar material. There is some evidence that the X-ray 
flux is mainly thermal, i ndicating the blas t-wave interac- 
ti on with dense matter (IPark et afll2007f) . 

iBerezhko fc K scnofontov (2006) compared radio and 
X-ray data with their calculations and concluded that 
there is a high (amplified) magnetic field inside the super- 
nova remnant. The predicted TeV gamma-ray flux, re- 
sulting mostly from the decay of neutral pions produced 
in interactions of the accelerated cosmic rays, is almost 
one order below the level obtained in previous searches 
by Imaging Atmospheric Cherenk ov Telescopes (IACTs) 
such as CANGARO O-II in 2001 ([Enomoto et al.l l200l 
and H.E.S.S. in 2003 (|Rowell et al]|2005fh The wide field 
of view searched by H.E.S.S. did, however, yield marginal 
excesses in a region north-east of SN 1987A and at LMC 
X-l. 

CANGAROO-III is one of two stereoscopic IACTs lo- 
cated in the southern hemisphere. With a three-fold co- 
incidence, the sensitivity is significantly improved from 
the previous single telescope (CANGAROO-II). We also 
increased the observation period by a factor of two in 
the 2006 observations compared to the 2004 and 2001 
seasons. We here report the results of these observations. 

The previous H.E.S.S. analysis considered not only SN 
1987A but also four other high energy objects within the 
same field of view ([Rowell et al.ll2005D . We follow that 
idea and present upper limits on the integral gamma-ray 
fluxes from the super-bubble 30 Dor C, the Crab-like pul- 
sar PSR B0540-69, the X-ray binary LMC X-l, and the 
supernova remnant N157B. Super-bubbles are expected 
to be efficient accelerators of cosmic rays (jParizot et al.l 
l2004h and a shell feature near the star 30 Dor C, within 
5 arcmin of SN 1987A, sho ws an indication o f non- 
thermal hard X-ray emission (jBamba et all [2004V The 
energy required to produce the bubble is estimated to be 
7xl0 51 erg, which is higher t han generally expected for 
a single supernova remnant (|Ueno. Bamba. fc Kovamal 
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I2003T) . The other three sources belong to categories 
of high-energy astrophysical objects that may acceler- 
ate particles to TeV energies, but, considering their 
large distance of 50kpc, they are less likely to be de- 
tected by current IACTs than their Galactic counter- 
parts. However, as they are well within the field of 
view of the CANGAROO-III telescopes, and given the 
marginal H.E.S.S. excess toward LMC X-l, it is straight- 
forward to search for any evidence of gamma-ray emis- 
sion from them. 

2. CANGAROO-III STEREOSCOPIC SYSTEM 

The CANGAROO-III stereoscopic system consists of 
four imaging atmospheric Cherenkov telescopes located 
near Woomera, South Australia (31°S, 137°E). Each 
telescope has a 10 m diameter segmented reflector con- 
sisting of 114 spheric al mirrors each of 80 cm diameter 
(|Kawachi et al.ll2001l ). providing a total light collecting 
area of 57.3 m 2 . The spherical segments are mounted on 
a parabolic frame with a focal length of 8 m. The first 
telescope, Tl, which w as the CANGAROO-II telescope 
(|Enomoto et al.ll2002al h was not used for these observa- 
tions due to its smaller field of view and higher energy 
threshold. The second, third, and fourth telescopes (T2, 
T3, and T4) were used for the observations described 
here. The camera syste ms for T2, T3, a nd T4 are identi- 
cal and are described in iKabuki et al . ( 2003). The tele- 
scopes are located at the east (Tl), west (T2), south 
(T3) an d north (T4) corners of a diamond with sides of 
~100m (jEnomoto et al.ll2002bf ). 

3. OBSERVATIONS 

The observations were carried out in the periods from 
2004 Nov 11 to 14 (MJD 53320-53323) and from 2006 
Dec 12 to 27 (MJD 54081-54096) using "wobble mode" 
in which the pointing position of each telescope was 
shifted i n declination betw een ±0.5 degree every 20 
minutes (jDaum et al.lll997T ) from the target: (RA, dec 
[J2000]) = (83.866139°, -69.269577°). We took no OFF 
source runs. The most sensitive region of the field of 
view is within one degree of the average pointing po- 
sition. LMC X-l is located 0.6 degrees south-east from 
SN 1987A. Considering our angular resolution of 0.24 de- 
gree, there is some overlap between SN 1987A, 30 Dor C, 
PSR B0540-69, and N157B. 

In the 2004 observation, data with GPS time stamps 
were recorded for T2, T3 and T4 individually when more 
than four photomultiplier (PMT) signals exceeded 7.6 
photoelectrons (p.e.). In the offline analysis stage we 
combined the data when the three telescope's GPS times 
coincided. In the 2006 observation, a hardware coinci- 
dence was used to se lect any two triggered telescopes 
(|Nishiiima et al.l |2005[ ) . The images in all three tele- 
scopes were required to have clusters of at least five adja- 
cent pixels exceeding a 5 p.e. threshold (offline three- fold 
coincidence). The event rate was reduced to 5^8 Hz by 
this criterion depending on the elevation angle and year. 
Based on time dependence of these rates we can remove 
data taken in cloudy conditions. The effective obser- 
vation times for 2004 and 2006 were 632 and 1316 min, 
respectively. The corresponding mean zenith angles were 
42.1° and 40.0°. 

The light collecting efficiencies, including the reflec- 
tivity of the segmented mirrors, the light guides, and 



the quantum efficiencies of the photo multiplier tubes 
were m onitored by a muon-ring analysis (jEnomoto et al.l 
l2006al) with the individual trigger data in the same peri- 
ods. The light yield per unit arc-length is approximately 
proportional to the light collecting efficiencies. Deteri- 
oration is mostly due to dirt and dust settling on the 
mirrors. 

4. ANALYSIS 

Here, we briefly describe the analysi s procedures, 
which are identical with those described inlKabuki et al. 
(|2007l). Mo re details can be foun d in lEnomoto et al. 
(j2006;Ji and lEnomoto et al.l (l2006bf). 

First, the image moments (Hillas 1985) were calculated 
for the three telescopes. The incident gamma-ray direc- 
tion was determined by minimizing the sum of squared 
widths (xq- weighted by the photon yield) of the three 
images seen from the assumed position (fitting parame- 
ter) with a constraint on the distances from the intersec- 
tion point to each image center. 

In order to derive the gamma-ray likeli ness, we used 
the Fisher Discrimina nt (hereafter FD) ((Fisher 1936; 
lEnomoto et al.l l2006a). Input parameters were 

P = (W2, W3, WA, L2, Li, LA), 

where W2, W3, W4, L2, L3, LA are energy-corrected 
widths and lengths for the T2, T3, and T4 and assume 
that a linear combination of 

FD = a ■ P, 

provides the best separation between signal and back- 
ground, then the set of linear coefficients (a) should be 
uniquely determined as 

rfmg - Msg 
a = — — , 

E s ig + Ebg 

where fl is a vector of the mean value of P for each sam- 
ple and E is their correlati on matrix. We previou sly used 
it in Vela Pulsar analysis (|Enomoto et al.|[2006af ) to sep- 
arate "sharp" (gamma-ray-like) images from "smeared" 
(background) ones. The values of fl s i g , Peg? E sig , and 
Ebg can be calculated from the Monte-Carlo and obser- 
vational data (OFF-source runs), respectively. 

We rejected events with any hits in the outermost layer 
of a camera ("edge cut"). These rejected events suffer 
from finite deformations, especially in the length distri- 
bution, which results in deformations of the FD. 

We then derived FD distributions position by posi- 
tion. Comparing those in the signal region and control 
background region, we can determine the gamma-ray- 
like events. Here, we use Monte-Carlo simulations to de- 
termine the FD distribution of gamma-ray signals. Note 
that in the gamma-ray simulations we used a spectrum 
proportional to E~ 21 . The fit of the FD distribution 
of source position with the above emulated signal and 
control background functions were carried out, to derive 
the number of gamma-ray-like events. This is a one- 
parameter fitting with the constraint that sum of signal 
and background events corresponds to the total number 
of events. These coefficients can be derived exactly ana- 
lytically. 
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5. RESULTS 

The threshold of this analysis is considered to be 
~1 TeV which is higher than the typical CANGAROO- 
III threshold due to the larger zenith angle of the observa- 
tions, of around 40 degrees. In order to derive morphol- 
ogy, we segmented the field of view into 0.2 x 0.2 degree 2 
square bins. The FD distributions for corresponding 
bins are made and fitted. The control-background re- 
gion is defined as the second closest layer, which is more 
than 0.3 degree from the center of target region, i.e., 
larger than the 0.24 degree point spread function (PSF). 
The statistics of the control-background is, therefore, six- 
teen times larger than that of signal bin. The results 
are shown in Figs. [TJ The left panel is obtained from 
the 2004 data and the right from 2006. The average 
telescope pointing position is indicated by the cross cen- 
tered and labeled as "SN 1987A", with the other four 
high energy objects within this field of view indicated 
by crosses. Our sensitivity is considered to be limited 
up to one degree in radius from the center. The PSF is 
considered to be a 0.24 degree radius circle. The signif- 
icance distributions (excess divided by the statistical er- 
ror) are approximately normal (Gaussian) distributions. 
The best fit Gaussians have mean values of — 0.31±0.08 
(2004) and -0.18±0.09 (2004), and standard deviations 
of 1.21±0.06 (2004) and 1.31±0.06 (2006), within sys- 
tematic uncertainties. 

At first we investigate the region of SN 1987A. Note 
that this region at some energies is contaminated by 
30 Dor C, a highly energetic source. To indicate the ex- 
cess distribution from its center, we show so-called B 2 dis- 
tributions in Figs. [2] Here, the control-background sam- 
ple was selected in the region 6* 2 =(0. 1-0. 2) [degree 2 ]. The 
background-subtracted signals are shown. Both 2004 and 
2006 data have statistically insignificant excesses near 
8 2 = 0. The x 2 /DOFs (degree of freedom) for null as- 
sumptions are 18.6/15 and 13.2/15 for 2004 and 2006, 
respectively. These statistically insignificant excesses ap- 
pear as the small excess count peaks around the centers 
of the plots in Fig. [TJ Note that our PSF corresponds to 
6 2 <0.06 degree 2 . We, therefore, proceed assuming that 
there is no signal present. In order to show the consis- 
tency between 2004 and 2006, the vertical unit is unified 
to count rate. 

In order to determine whether or not there is a gamma- 
ray excess around SN 1987A, we made the FD distribu- 
tions within the PSF. The control-background was made 
from the region 6* 2 =(0.1-0.2). The fitting results are 
shown in Fig. [3] The background level differs between 
the two years. This is due to sky conditions, mirror re- 
flectivity, time-dependence of the blur spot size (which 
improves after maintenance periods), etc. These are con- 
sidered to be the major systematics for estimating the 
observational gamma-ray flux in this work. For the light 
correction efficiency, we may have 10% uncertainty at 
maximum. 

Again there is no statistically significant excess. We 
carried out a similar analysis, year by year, energy 
threshold by energy threshold. The obtained 2a upper 
limits (ULs) are listed in Table [TJ Here we used a E~ 21 
spectrum for the gamma-ray simulation. The ULs range 
from 6-17% crab. The worse upper limits in the lower 
energy range originated from the statistically insignifi- 



TABLE 1 

The 2a upper limits to the integral fluxes from SN1987A 

AT TWO EPOCHS WITH THREE DIFFERENT ENERGY THRESHOLDS. 



Year Excess Upper Limit Energy Threshold Flux Upper Limit 
Events GeV cm 2 s _1 



2004 35 1200 1.7 x lO -11 

2004 21 2870 5.7 x 10" 13 

2004 4.6 7180 1.0 x 10" 13 

2006 116 1080 2.7 x 10~ 12 

2006 19 2550 2.5 x 10" 13 

2006 5.3 6350 6.1 x 10" 14 



cant excess around the center of the field of view. At 
higher energies, we do not see any excess. 

The four other high energy sources in the surrounding 
region, 30 Dor C, PSR B0540-69, LMC X-l, and N157B, 
can be analyzed in the same way. They are all within 
a one degree circle from the average pointing position 
(the sensitive region) . The summary is listed in Table [H 
These are two year averages. Note that, considering our 

TABLE 2 

The 2cr upper limits to the integral fluxes at energy 

GREATER THAN 1120 GeV FOR FIVE HIGH ENERGY OBJECTS. THE 
2004 AND 2006 DATA WERE AVERAGED. 



Target 


Excess Upper Limit 


Flux Upper Limit 




Events 


— 9 — 1 

cm z s 1 


SN 1987A 


126 


1.8 X 10" 12 


30 Dor C 


165 


2.5 x 10" 12 


PSR B0540-69 


66 


1.1 x 10~ 12 


LMC X-l 


74 


1.1 x 10~ 12 


N157B 


167 


2.4 x lO" 12 



PSF of 0.24 degree, SN 1987A and 30 Dor C are highly 
confused, with other pairs being less confused. Also some 
of the control-background regions are confused. The ac- 
ceptance is a slowly and smoothly decreasing function 
and is reduced to 65% at a distance of one degree from 
the center. 

6. DISCUSSION 

Our morphology indicates a negligibly small excess 
near the center of field of view both in the 2004 and 
2006 data. This is statistically <2a each year. We note 
that a similar feature can also be seen in the H.E .S.S. 
data shown in Figs. 2 and 3 of iRowell et alj (|2005h . A 
deep survey, i.e., longer observation, around this region 
is awaited. The region around LMC X-l is consistent 
with zero in our data. 

The upper limits of SN 1987A in Table[TJare plotted in 
Fig. [H The results presented here (in blu e and black) are 
comp ared with CANGARO O-II (red! (lEnomoto et all 
f2003h and H.E.S.S. (green) (|Rowell et al] 12005ft result"^ 
The upper limits in this work are improved by fac- 
tors of up to 10 on the 2001 result. That is con- 
sistent with the longer observation time and better 
signal-to-noise ratio resulting from the three-fold co- 
incidence. However, the theoretical prediction even 
7300 days after the superno va (the dashed curve) 
(|Berezhko fc Kseno fontov 2006) is still a factor of three 
below th e results of the observations. No te that the pre- 
dictions (|Berezhko &: Ksenofontov! I2006f ) were made in 
the frame of spherical symmetric approach. The exis- 
tence of strongly asymmetric interstellar matter struc- 
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PSR B0540-69 
0. 
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Fig. 1. — Excess count maps. The left panel is obtained from the 2004 data and the right from the 2006 data. The average pointing 
center is SN 1987A itself. The positions of four other high energy objects are also indicated by the black c rosses with labels. The PSF 
is shown by the black d ashed circle. The white contours are derived from 4.85 GHz PMN survey radio data {Wright et al. 1994) obtained 
from Skyview (NASA) UNASA II2007T ). 




S 0.4 



(degree ) 



Fig. 2. — Distributions of 8 2 [degree 2 ]. The left panel is obtained 
from the 2004 data and the right from the 2006 data. The vertical 
scale is normalized to counts/minute per bin [0.2 degree 2 ]. The 
background-subtracted signals (obtained by the fitting procedure 
described in the text) are plotted. 

ture, which is a dense inner ring, makes their prediction 
of gamma-ray flux uncertain. According to their rough 
estimation, this uncertainty is not very large (about a 
factor of 2) due to the stronger supernova shock decel- 
eration in the denser medium. The solid curve is their 
prediction for 8249 days after the supernova, still lower 
than the upper limits obtained in this work. In order to 
constrain this model, we need further improvements on 
the hardware and/or analysis methods. A deeper obser- 
vation by H.E.S.S., such as T>30 h, around -2010 might 
give conclusive results. Also proposed future large-scale 
IACT arrays such as CTA ([Hermann 1 |2007| ) would be 
very promising. 

The limit obtained for 30 Dor C in Table [5] is about 
15% crab level at —1 TeV. Considering the distance of 
50kpc to LMC, this upper limit corresponds to a huge 
cosmic-ray density comparing with that of the Crab neb- 
ula (distance of 2kpc). Even if there is supernova com- 




FlG. 3. — FD distributions. The left panel is obtained from 
the 2004 data and the right the 2006 data. The vertical scale 
is counts/minute per bin [arbitrary unit]. The black crosses are 
obtained from the on-source region (centered around SN 1987A 
within the PSF), the green histograms the best fit of the control 
backgrounds, the blue crosses the background-subtracted signals, 
and the red histograms the best fit signal distributions. 

plex there, the expected TeV gamma-ray flux is still be- 
low our sensitivity limit. Increasing the cosmic-ray target 
density such as either inter-stellar-matter or photon den- 
sity within plausible ranges, the expected flux does not 
exceed this high upper limit. According to the standard 
es timation on the hadronic gamma- r ay production such 
as iDrurv. Aharonian. fc Volki (|1994D : iNaito fc Takaharal 
(|1994h with an assumption of the target density of 
lp/cm -3 , a total kinetic explosion energy of 7 x 10 51 ergs 
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Fig. 4. — Spectral energy distribution. The blue points and line 
are obtained from the 2004 data and the black from the 2006 data. 
The r ed points are the CA NGAROO-II upper limits from 2001 
data {Enomoto ct al. 2003). The green points are the H.E.S.S. 
data from 2003 (Rowcll et al.||2005|) . The lines are obtained from 
Fig. 4 of Bcrczhko & Kscnofontov (2006). The solid and dashed 
curves are the predictions of the gamma-ray flux 8249 and 7300 
days after the supernova, respectively. 

gives a 1 milli-crab flux of gamma-rays at the solar sys- 
tem. Again future larger-scale projects are required to 
probe such flux levels. 

For the other three high energy objects, we set upper 
limits far below those of previous observations. We, how- 
ever, could not find any physically important constraints 
on their activity at high energies again due to the dis- 
tance of 50kpc. 



7. CONCLUSION 

We have observed SN 1987A at two epochs: MJD 
53320-53323 (2004 Nov) and MJD 54081-54096 (2006 
Dec), approximately 6500 and 7200 days after the super- 
nova, respectively. The effective observation times are 
10.5 and 20 h, respectively. No statistically significant 
gamma-ray signals were obtained and we set 2er upper 
limits on the integral fluxes of 6-17% crab at ~1 TcV. 
These are slightly lower than the previous CANGAROO- 
II and H.E.S.S. upper limits. Although theory predicts a 
factor of >3 lower flux 7300 days after supernova, it also 
predicts that the future flux level might exceed the sensi- 
tivity limit of the existing and future arrays. Continued 
monitoring of this object at TeV energies is still therefore 
meaningful in constraining theoretical models. It is also 
important to improve the sensitivity of the o bservation 
meth ods, as proposed by the CTA project ((Hermann] 
[2001 . 

In addition, we set 2a upper limits for four more high 
energy objects inside our field of view: 30 Dor C (super 
bubble), PSR B0540-69 (Crab-like pulsar), LMC X-l 
(X-ray binary), and N157B (supernova remnant). 
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